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Abstract. We report on observations of correlated behavior between the prompt y-ray and optical emission from GRB 
0803 19B, which (i) strongly suggest that they occurred within the same astrophysical source region and (ii) indicate that 
their respective radiation mechanisms were most likely dynamically coupled. Our preliminary results, based upon a new 
cross-correlation function (CCF) methodology for determining the time-resolved spectral lag, are summarized as follows. 
First, the evolution in the arrival offset of prompt y-ray photon counts between Swift-BAT 15-25 keV and 50-100 keV energy 
bands (intrinsic y-ray spectral lag) appears to be anti-correlated with the arrival offset between prompt 15-350 keV y-rays 
and the optical emission observed by TORTORA (extrinsic optical/y-ray spectral lag), thus effectively partitioning the burst 
into two main episodes at --^ T + 28 ± 2 sec. Second, prompt optical emission is nested within intervals of (a) trivial intrinsic 
y-ray spectral lag (~ T + 12 ± 2 and ~ T + 50 ± 2 sec) with (b) discontinuities in the hard to soft evolution of the photon index 
for a power law fit to 15-150 keV Swift-BAT data T + 8 ± 2 and ~ T + 48 ± 1 sec), both of which coincide with the rise 
(~ r + 10 ± 1 sec) and decline (~ T + 50 ± 1 sec) of prompt optical emission. This potential discovery, robust across heuristic 
permutations of BAT energy channels and varying temporal bin resolution, provides the first observational evidence for an 
implicit connection between spectral lag and the dynamics of shocks in the context of canonical fireball phenomenology. 
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1. INTRODUCTION 

Swift's unique dynamic response and spatial localization precision, in conjunction with correlative ground-based 
follow-up efforts, has resulted in the collaborative broad-band observations of GRB 0803 19B [1]. In the context of 
an analysis focused on confronting the lag-luminosity relation [2] in the Swift era [3], a correlation was observed 
between the evolution of time-resolved spectral lag and the behavior of the extraordinarily well-sampled prompt 
optical emission light curve associated with GRB 0803 19B. In general, the spectral lag is determined via either a 
peak pulse fit |4] or cross-correlation function (CCF) analysis |5]. Previous studies have reported on the variability of 
spectral lag throughout burst emission |@], as well as its correlation to pulse evolution |7]. In this work, we develop 
a new method to calculate the time-resolved spectral lag via a modification to the traditional CCF approach. In this 
manner, we are able to explore for the first time the evolution and correlation of the intrinsic y-ray spectral lag with 
prompt optical emission. We interpret these correlated behaviors as strong observational evidence that the prompt 
optical and y-ray emission of GRB 0803 19B took place within the same astrophysical source region, with indications 
that their respective radiation mechanisms where dynamically coupled throughout the prompt phase of the burst. 



' Adapted from a contribution to the Proceedings of the 2008 Nanjing GRB Conference. Edited by Y. F. Huang, Z. G. Dai and B. Zhang. 
^ Correspondence toMichael.Stamatikos-l@nasa.gov. 



2. METHODOLOGY 



The Burst Advocate (BA) script was used to generate BAT event data and quality maps, which were used to construct 
32 ms binnecQ (mask-weightecQ background-subtracted) Ught curves (temporal photon spectra) for Ea = 15 — 25 keV 
and Eg = 50— 100 keV energy bandfl via the BATBINEVT analysis task. The resultant FITs files were analyzed 
via the CCF method in order to quantify the temporal correlation between the two series of GRB light curves in 
differing energy channels, i.e. Ea and Eb, via the following: 

where ffy^,, = (^|„)'''^- The CCF between Ea and Eb peaks at a given temporal offset (fo), known as the spectral 
lag (tab), which is defined as positive if the systematic shift in the arrival times of photon counts between pairs of 
light curves results in higher energy {Eb) photons arriving before those of lower energy [Ea). The generic CCF, based 
upon the Pearson Correlation, is defined in Equation[T] If the mean is subtractec^, then v(f) — d{t) — {d{t)). If the 
mean is not subtracted, then v{t) — d{t) — b{t), where an attempt is made to remove background counts b{t). Our tests 
indicated that both methods agreed when one calculates the time-averaged spectral lag, i.e. the spectral lag determined 
over the entire duration of the GRB. This included a test of subtracting the mean of a light curve that was nested within 
two large background intervals, which had the effect of reducing the mean and increasing the correlation amplitude. 
Note that if one takes an infinite interval, the mean would go to background, which fluctuates about zero, resulting in 
effectively not subtracting the mean in the first place, as prescribed for transient sources, i.e. GRBs [5]. 

However, when one tries to extract a time- resolved spectral lag, i.e. the spectral lag over a segment of the light curve, 
then only the mean subtraction method consistently works. This is a consequence of edge effects introduced by the 
assumptions of a stationary versus a transient temporal signal characterization. If one subtracts the mean, there is an 
assumption that the signal exists as a stationary wave throughout time, i.e. beyond the sampled time series. If one does 
not subtract the mean, then the assumption is that the signal is entirely contained within the time series. Both methods 
disregard data beyond the overlapped interval due to the introduction of a spectral lag shift. 

Hence, for a bright GRB like 0803 19B, a non-mean subtracted time-resolved spectral lag analysis is dominated by 
the segment of the light curve described by a square wave, and the real light curve signal is treated as fluctuations. A 
mean-subtracted time-resolved spectral lag analysis will extract the spectral lag since it smoothes out the amplitude 
and focuses on the signal variance. Hence, we utilized the mean-subtracted definition of the CCF, since it is both more 
conservative and robust. Although it has been demonstrated that spectral lag variability is ubiquitous in GRBs 
previous methods required relatively well-behaved (FRED-like) sub pulses. Our new methodology treats structure 
without regard to its functional form and affords greater flexibility for time-resolved spectral lag evolution studies. 

In order to test the viability of an interval, a large known spectral lag was artificially introduced, which shifted one 
of the test light curve pairs. Intervals where the the analysis was unable to recover the artificial spectral lag due to 
either low SNR or lack of structure were not used. The CCF method, regardless of definition, is sensitive to the SNR 
and, more importantly, to the shape of the region. Morphological tests using simple square, triangular and Gaussian 
pulse shapes have illustrated that both mean and non-mean subtracted CCF definitions failed to recover the artificial 
spectral lag that was obscured in regions where (i) only slopes were compared, (ii) SNR was low and/or (iii) error bars 
were large. This puts a fundamental limit on determining the time-resolved spectral lag for discrete intervals, which 
precluded an arbitrary systematic sampling of GRB 0803 19B, as illustrated in Figure [T] Panel D. Hence, although 
our methodology is robust, it does require a selection on intervals with structured variability. This is not an issue for 
time-averaged results, since the burst's structured variability interval is nested between periods of quiescence. 



' A binning of 32 ms was used based upon a heuristic exploration of various bin resolutions (e.g. 4 ms, 8 ms, 16 ms, 32 ms, 64 ms, etc.). Once a 
spectral lag measurement was made, it was reproduced using light curves of several adjacent bin resolutions. In principle, the binning of the light 
curve should be smaller than the spectral lag in order to resolve the CCF peak. 

Although one increases the signal-to-noise ratio (SNR) using non-mask-tagged (raw, non-background subtracted) hght curves, artifacts either 
intrinsic to the detector or associated with the background may mimic signal pulses and thus be treated as such in a CCF analysis. Also, the 
background may not be trivial to model since it may be highly variable due to slewing. Hence, in order to be conservative and consistent with future 
studies, background subtracted, i.e. mask-tagged, light curves were used throughout the analysis. 

' Permutations of canonical BAT channel (1-4) light curve pairings were investigated. Ultimately BAT channels 1 (15-25 keV) and 3 (50-100 keV) 
were used in the analysis since they represented the largest energy differential with the highest count rate significance. 

* Mean subtraction is the formalism identical to that used by the C_CORRELATE function of IDL 6.2. The formal definition may be found onhne 

(http : //cl . dmf . arm. gov/base/idl_6 . 2 /C_CORRELATE . html). 



The error bars on the CCF function were determined via a primary Monte Carlo (MC) simulation, which were 
constructed from 1000 realizations of light curve pairs, within the la error bars of the data. The distribution of 
CCF values within a given bin defined its error bar Once the error bars were fixed in this manner, a secondary MC 
simulation constructed a subsequent set of 1000 realizations of light curve pairs, within the la error bars of the data. 
Each resulting CCF (with errors fixed from the primary MC simulation) was subjected to a Gaussian fit: 

3;=Ae-("^)'/(2'^'), (2) 

whose peak was extracted via the IDL function MPFITPEAjfl. The peak fit to the data was used to determine the 
spectral lag value, while the numerical standard deviation of the distribution of simulated peaks was used to determine 
the error. We interpret this bootstrapped result as the la spectral lag error, as illustrated in Figure[T] Panel D. 

In general, spectral lag calculations are sensitive to a series of selection parameters that include the energy band 
pass of each comparative light curve, temporal bin resolution and emission interval. The CCF and pulse fit methods 
are both based upon a functional expression for the peak fit (Gaussian, cubic, etc.) while the former requires a bin shift 
range. A self-consistent comparison, which includes uncertainties, must account for these selection effects. Efforts 
to standardize the spectral lag analysis are ongoing and would facilitate direct comparative studies across detector 
sensitivities and energy band passes. 



3. RESULTS & DISCUSSION 

The generic agreement of the overall temporal coincidence and morphology between the prompt 7-ray and optical 
light curves (Figure [T] Panels A & B) suggests that they arose from a common source region. However, separate 
radiation mechanisms were most likely responsible since the extrapolated y-ray flux density to the optical band was 
deficient by ~4 orders of magnitude when compared to observation y], Igl |9|] . The steep rise/decline, short duration 
and lack of increasing pulse width of the prompt optical emission disfavor external forward/reverse shocks. Hence, 
internal shocks have been suggested as the source region with synchrotron emission responsible for the optical and 
inverse Compton scattering/synchrotron self Compton for the 7-rays, with associated GeV photon emission [L .8]. 
Alternatively, it has been suggested that non-relativistic forward internal shocks generated the prompt optical emission, 
while relativistic reverse internal shocks generated the prompt 7-rays, with sub-GeV/MeV photon emission 1911. S uch 
high energy emission may be tested by Fermi (formally known as GLAST) via joint analyses with Swift-BAT flOtl . 

Our preliminary results are illustrated in Figure[T] Panels A-D. There are several observations that lead to effectively 
separating the burst's duration into two main episodes partitioned roughly at the midpoint of ~ J -|- 28 ±2 sec. The 
first is that the bimodal evolution of the intrinsic 7-ray spectral lag increases at t >^ T + 28 ± 2 sec, which appeal's 
to be anti-correlated with the extrinsic (optical/7-ray) spectral lag observed via the smoothed Gaussian fits, where 
the spectral lag is of the order of a few seconds when t<~r+28±2 sec (Panels A & B). Beyond this common 
midpoint, the optical and 7-rays do not correlate as well or at least are ambiguously correlated, i.e. extrinsic (optical/7- 
ray) spectral lag is either zero or negative at later times (ambiguity due to peak misalignment). Hence, the intrinsic 
time-resolved 7-ray spectral lag is maximum att>^r + 28±2 sec, while the extrinsic time-resolved (optical/ 7-ray) 
spectral lag is maximum att<^r + 28±2 sec. In addition, an independent analysis lITlll of BAT 15-150 keV light 
curves has revealed that the characteristic variability timescale of GRB 080319B was ^100 ms for t < ^ rH-28 sec 
and '^l sec for t > ^ rH-28 sec. Furthermore, time-resolved spectral analysis by Konus-Wind illustrated that Epg^k 
decreased from 75 1 ± 26 keV to 537 ± 28 keV at - T + 24 ± 2 sec lH . 

The hard to soft evolution of the photon index for time-resolved power law fits to 15-150 keV Swift-BAT data occurs 
in steps at ^ T + 8 ± 2 and ^ T + 48 ± 1 sec, which coincide with the respective rise (^ r + 10±l sec) and decline 
('^ r + 50 ± 1 sec) of prompt optical emission (see Panels A & C). This is consistent with Konus-Wind time-resolved 
spectral analysis and hardness ratios |1]. We also note that intervals of trivial intrinsic spectral lag coincide with the 
rise and decline of prompt optical emission (see Panels A & D), at~r+12±2 and ^ T + 50 ± 2 sec, respectively. We 
interpret these correlated behaviors as strong observational evidence that the prompt optical and 7-ray emission took 
place within the same astrophysical source region, which until now has only been conjecture. We also find indications 
for a dynamical coupling between the radiation mechanisms, perhaps via the processes mentioned above llIlH,^]- 



' See http :/ / astro . berkeley . edu/~heiles /idlprocs /mpf it /mpf itpeak . pro, for more information. 




Seconds Since Swift-BAT GRB 080319B Trigger Time (06:12:49 UT on IVlarch 19, 2008) 



FIGURE 1. Panel A - Prompt optical flux density as observed by TORTORA for ~1.3 second exposure intervals (black circles). 
A cumulative fit (solid red line) is based upon 3 Gaussian dashed curves with peaks at ~T+19 sec (black), ~T+34 sec (blue) and 
~T+45 sec (magenta). Data has been digitized from Figure 1 of U^. Panel B - Swift-BAT prompt y-ray (32 ms bin) light curves 
for canonical energy channels 1 (15-25 keV, solid black) and 3 (50-100 keV, solid blue). 15-350 keV summed, 1.3 second bin light 
curve illustrated in solid green. A cumulative fit (solid red line) is based upon 3 Gaussian dashed curves with peaks at ~T-l-17 sec 
(black), ~T-l-40 sec (blue) and ~T-l-49 sec (magenta). Panel C - Photon index for power law fit to 15-150 keV Swift-BAT 

data, with 1 o error bars. Note that hard to soft evolution occurs in discontinuous steps at ~ 7 + 8 ± 2 and ~ J -I- 48 ± 1 sec, which 
coincide with the respective rise (~ T + 10± 1 sec) and decline (~ r + 50± 1 sec) of prompt optical emission (see Panel A). Panel 
D - Time-resolved, intrinsic spectral lag between BAT 32 ms bin light curves for 15-25 keV and 50-100 keV energy bands, based 
upon a Gaussian fit to the CCF peak, with IcT error bars derived from 1000 Monte Carlo bootstrapped simulations. Intervals of 
trivial spectral lag occur at ~ T + 12 ± 2 and ~ T + 50 ± 2 sec and also coincide with the rise and decline of optical emission. 
Overall, the temporal evolution of the time-resolved spectral lag can be described by two main episodes, where the intrinsic (BAT 
Channel 1/Channel 3) y-ray spectral lag is greater for t>~r + 28±2 sec. Note that this is anti-correlated with the extrinsic 
optical/y-ray spectral lag observed via the smoothed Gaussian fits to the optical emission, where the spectral lag is of the order of 
a few seconds att<~r + 28±2 seconds, and either zero or negative at later times (ambiguity due to peak misalignment). Such 
correlated spectral and temporal behavior indicates that the prompt optical and y-ray emission took place within the same source 
region, with radiation mechanisms that were most likely dynamically coupled. 



This potential discovery, provides the first observational evidence for an implicit connection between spectral lag 
and the dynamics of internal shocks in the context of canonical fireball phenomenology. A full theoretical analysis 
of this result is currently in preparation. Future work includes an application of our methodology to observations of 
a subset of bursts with prompt optical emission to probe either the unique or ubiquitous nature of these observations. 
Ultimately, understanding the mechanism(s) responsible for spectral lag may reveal a fundamental and unprecedented 
view from within the GRB fireball. 
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